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The effective on-site Coulomb interaction (Hubbard U) between localized d electrons in 3(i, 
Ad, and 5rf transition metals is calculated employing a new parameter- free realization of the con- 
strained random-phase approximation using Wannier functions within the full-potential linearized 
augmented-plane-wave method. The U values lie between 1.5 and 5.7 eV and depend on the crystal 
structure, spin polarization, d electron number, and d orbital filling. On the basis of the calcu- 
lated U parameters, we discuss the strength of the electronic correlations and the instability of the 
paramagnetic state towards the ferromagnetic one for 3d metals. 
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Density functional theory (DFT) within the local- 
density approximation (LDA)i is a reliable method for 
calculating ground-state properties of solids with weak 
electronic correlations, i.e., for U /W < 1, where U is the 
effective on-site Coulomb interaction between localized 
electrons (Hubbard U) and W is the bandwidth. How- 
ever, LDA often fails to describe systems with intermedi- 
ate {U/W ^ 1) and strong {U/W > 1) electron correla- 
tions, such as transition-metal oxides, rare earths, Kondo 
systems, etc. Properties of these materials are usually 
calculated with phenomenological many-body Hamilto- 
nians such as the Hubbard^ or the Anderson impurity 
model.™ However, in these models the Coulomb and also 
the one-particle hopping matrix elements are typically 
empirical parameters that are determined such that the 
employed model reproduces experimental results of in- 
terest. 

For a long time, DFT-LDA and many-body model 
Hamiltonian methods have been separate and comple- 
mentary approaches. This has drastically changed with 
the advent of the dynamical mean-field theory (DMFT)i 
that merged with LDA to a novel computational method 
referred to as LDA-|-DMFT)^ which developed to a mod- 
ern many-body approach for treating correlated electron 
materials. In retrospect, the so-called LDA+U method,- 
an early attempt correcting the LDA functional by in- 
troducing a simple mean-field-like Hubbard U term for 
localized d or f states, and today routinely applied to a 
broad spectrum of systems can be regarded as its static 
limit. Both LDA+U and LDA-I-DMFT as well as other 
approaches not mentioned here have in common that 
they rely on the Hubbard U as an additional parame- 
ter. Frequently the exact value of U is unknown which 
impedes the predictive power of these approaches. 

The problem of calculating the Hubbard U parameter- 
free for transition metals (TMs), i.e., from first princi- 
ples, has been addressed by several authors i^""— A num- 
ber of different approaches have been proposed. Among 
them, constrained local-density approximation (cLDA)S 
is the most popular one. However, cLDA is known to 
give unreasonably large U values for late TMs due to dif- 



ficulties in compensating the self-screening error of the 
localized electronsJ^ Furthermore, the frequency depen- 
dence of U is unattainable. On the other hand, the con- 
strained random-phase approximation (cRPA), though 
numerically much more demanding, does not suffer from 
these difficulties. In contrast to cLDA, it also allows to 
access individual Coulomb matrix elements, e.g., on-site, 
off-site, intra-orbital, inter-orbital, and exchange. 

The aim of this Rapid Communication is to present 
a systematic study of the effective on-site Coulomb in- 
teraction (Hubbard U) between localized d electrons in 
TMs determined by means of first-principles calculations. 
Previous cRPA studies of U in TMs have focused only 
on the non-magnetic (NM) state of the 3d series and 
the results appeared to be strongly dependent on the pa- 
rameters used in the cRPA schemes J^iii In the present 
work we propose a new simple parameter-free cRPA ap- 
proach within the full-potential linearized augmented- 
plane-wave (FLAPW) method using maximally localized 
Wannier functions (MLWFs)J^ In addition to NM and 
magnetic states of the 3d series we consider dd and 5d 
elements in the periodic table. We discuss the strength 
of the electronic correlations and the instability of the 
paramagnetic state towards the ferromagnetic one for 3d 
TMs on the basis of the calculated Hubbard U parame- 
ters. 

The basic idea of the cRPA is to define an effective in- 
teraction U between the localized d electrons by restrict- 
ing the screening processes to those that are not explicitly 
treated in the effective model Hamiltonian.^*' To this end, 
we divide the full polarization matrix P = Pd + Pr, where 
Pd includes only d-d transitions and Pr is the remain- 
der. Then, the frequency-dependent effective Coulomb 
interaction is given schematically by the matrix equation 
U{u!) = [1 — vPr{uj)]~^v, where v is the bare Coulomb 
interaction and U{uj) is related to the fully screened in- 
teraction by U{uj) = [1 - C/(cj)Pd(w)]~iC/(cj). The static 
limit of the average diagonal matrix element of t/(u; — >■ 0) 
represented in a local basis can be regarded as the Hub- 
bard U parameter 

Although, cRPA is a general approach, its application 
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to materials with entangled bands is not straightforward. 
In these materials the localized d states that span the 
model subspace mix with extended s and p states, and 
there is no unique identification of the d-d transitions for 
constructing Pd- Several procedures have been proposed 
in the literature to overcome this problem. Aryasetiawan 
et al,^'^ suggested to use an energy window or a range of 
band indices to define the d subspace. However, the re- 
sults depended strongly on the chosen window or band in- 
dices. An alternative approach, in which the hybridiza- 
tion of the d states was switched off, was not burdened 
by additional parameters, but the U values turned out 
to be unphysically large for materials with strong sp-d 
mixing, e.g., early TMs, and the unphysical suppression 
of hybridization is unsatisfactory. 

In the present work we propose a new parameter-free 
procedure where Pd is directly constructed from the def- 
inition of the d subspace. The latter is spanned by a set 
of MLWFs <i,(r) = ^ Ek,„^ 7^R,™„(k)^£™(r), where 
iV is the number of k points, T^.^^{'k) is the unitary 
transformation matrix, ^^j^iic) are single-particle Bloch 
states of spin a and band index to, and R is the atomic 
position vector in the unit cell. We now define Pd as 
the polarization function that is generated by all tran- 
sitions that take place within the d subspace. To de- 
termine with what probability this applies to a given 
transition between extended Bloch eigenstates f^„-^ 
Vk+qm'j which may be mixtures of sp and d states, we 
multiply the probability that the electron resides in the d 
subspace before the transition p^^^ = „ I^^r .„„(k)P 
with the corresponding probability after the transition 
and obtain P^mPk+qm' ^^'^ probability for the transi- 
tion itself. Thus, Pd is constructed from summing over 
all transitions in the Lehmann representation multiplied 
with these probabilities. In this way, the resulting effec- 
tive interaction U{u!) only depends on the MLWFs that 
span the d subspace and is basically independent of the 
used electronic structure method. 

In a more formal approach, one can define Pd 
as the density correlation function Pd{rt,r't') — 
-i{'i'o\f[ndirt),hd{r't')]\^o) with the Kohn-Sham deter- 
minant ^'o, the time-ordering operator T, and the Heisen- 
berg density operator h{rt) = nd{rt)+nr{rt) decomposed 
according to the d subspace and the rest. Without time- 
dependent external fields, Pd only depends on the time 
difference t — t'. A Fourier transformation then yields the 
Lehmann representation described above. 

The ground-state calculations are carried out us- 
ing the FLAPW method as implemented in the 
FLEUR codei^ within the LDA exchange-correlation 
potential.ii The MLWFs are constructed with the 
Wcinnier90 codeJ^^ We include six bands per TM 
atom in the construction of the MLWFs, i.e., the 
five d bands and the itinerant s band. The ma- 
trix elements of the effective Coulomb potential U 
in the MLWF basis are given by U^^^^^.^^^^{uj) = 
//<:R(rXR(r)C/(r,r';a.)u;f:i,(r')<R(r') d^r d^r' 
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FIG. 1: (Color online) Average bare on-site direct Coulomb 
matrix elements between the d orbitals for TMs. In the inset 
we show the results for exchange Coulomb matrix elements. 



The effective Coulomb potential U (r, r'; ui) itself is calcu- 
lated within the cRPA implemented in the SPEX codeJ^ 
(for further technical details see Ref.fisl) . We define 
the average on-site diagonal (direct intra-orbital) and 
off-diagonal (exchange inter-orbital) matrix elements of 
the effective Coulomb potential as U — ^ J2n^ C^r 



and J = 



(d) 



a/3 

.m,n(m#„) f^Rm„;„m- The average off- 
diagonal (direct inter-orbital) Coulomb matrix elements 
are given by the relation U' = U ~ 2J. Although the 
matrix elements of the effective Coulomb potential are 
formally spin-dependent due to the spin dependence of 
the MLWFs, we find that this dependence is negligible in 
practice. 

We start with the discussion of the wnscreened (bare) 
Coulomb interaction in the TMs. Figure[T] shows the av- 
erage bare on-site direct {V) Coulomb matrix elements 
for the 3c?, Ad, and 5d TM series in the NM state. In 
the inset we show the results for exchange ( Jf,) Coulomb 
matrix elements. Note that among the 3d series Fe, Co, 
and Ni are ferromagnetic (FM) while Cr orders antifer- 
romagnetically. Also Mn is FM in the bcc structure with 
a = 2.91 A. For these elements we find that matrix el- 
ements of the bare Coulomb potential for magnetic and 
NM states are nearly identical. Within each series both V 
and Jb increase monotonically with the d electron num- 
ber. This can be explained by the fact that, as one moves 
from the left to the right within one row of the periodic 
table, the nuclear charge increases and causes the rf-wave 
functions to contract, which gives rise to the observed 
trend for V and Ji,. On the other hand, the localization 
of the d electrons decreases within one column of the pe- 
riodic table from 3d to 5d elements. As a consequence, 
V and Jb decreases in the same direction. This decrease 
is more pronounced for late transition metals. 

Efficient sp screening in TMs significantly reduces the 
bare Coulomb interaction V. Calculated Hubbard U and 
J parameters for the NM state of the TMs are presented 
in Fig. [21 Results for the magnetic states of the 3d el- 
ements are also included. For comparison, matrix ele- 
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TABLE I: Hubbard U for eg and t2g orbitals (in eV) for bcc 
V, Nb, and Ta and fee Ni, Pd, and Pt. 
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FIG. 2: (Color online) Left panels: U and U for the TM 
series. With open and filled diamonds we show U and U for 
the magnetic state of 3ii TMs. For comparison, experimental 
results from Refs.flSj (stars), [13 (pluses), and [2ll (crosses) as 
available are given. Right panels: The same for J and J. 



ments of the fully screened Coulomb interaction {U, J) 
are given. As seen in Fig. [21 in contrast to the bare direct 
Coulomb interaction V, the Hubbard U shows a non- 
monotonic behavior, i.e., it increases from the early TMs 
and reaches a plateau-like behavior around half-filling, 
whereas U is almost constant across the TM series, ex- 
cept for the elements with completely filled d shells like 
Cu. This behavior of the Hubbard U reflects a sub- 
stantial contribution of the d-d transitions to the fully 
screened Coulomb interaction U, especially around half- 
filling. In difference to the bare Coulomb V, the Hubbard 
U parameter is very sensitive to the d electron number 
and d orbital filing. In metals we are in the strong cou- 
pling limit, v\Pr\ ^ 1, and thus [/ ~ — p-. Since Pr 
depends mostly on the electronic structure of the screen- 
ing electrons, this explains why isovalent TMs with the 
same crystal structure of 4rf and 5d exhibit very similar 
U values but different values of Mn, Fe, and Co, which 
have different crystal structures. Furthermore, by con- 
strained NM and proper spin-polarized treatments of the 
magnetic elements, we show that spin polarization has a 
strong influence on U and U. The calculated U values for 
Cr, Fe, Co, and Ni turn out to be larger in the magnetic 
state than in the NM one, while the situation isjust the 
opposite for Mn, which was explained in Ref.iis! for the 
case of U by the different screening due to the available 
electrons at the Fermi energy. The same discussion holds 
also for partially screened U. Our calculated Hubbard U 
parameters for the 3 c? series are in good agreement with 
recent cRPA studies of Miyake et ali^ for late TMs as 
well as cLDA calculations of Nakamura et al^ for early 
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TMs. Experimentally, the Hubbard U parameters for 3d 
TMs are deduced from a combined use of Auger and x- 
ray photoemission spectroscopy. ^^"^^ Results from three 
different groups are included in Fig. [5] for comparison. 
As seen, the experimental U parameters are rather scat- 
tered. Our calculated U values are in good agreement 
with measurements of Kaurila et a/.— as well as Yin et 
al.^'^ for late TMs. So far, we have focused only on the 
effective intra-orbital direct Coulomb interaction U. The 
same discussion holds also for the inter-orbital direct and 
exchange Coulomb interaction U' and J, respectively. 
Note that in contrast to U and U' , renormalization of 
the J is rather small, i.e., J is close to the atomic value 
Jb- However, the d-d transitions substantially reduce J, 
especially for late TMs. It should also be noted that 
while the bare V has a long-range behavior, the U shows 
much faster damping. The calculated nearest-neighbor 
U values lie between 0.1 and 0.4 eV being maximal for 
TMs with half-filled d bands. 

In Table[T| we present the orbital dependence of the 
Hubbard U for bcc V, Nb, and Ta and fee Ni, Pd, and Pt. 
Fig. [31 shows the frequency dependence of U and J for the 
same elements. As can be seen, the crystal structure has 
a sizable influence not only on the orbital anisotropy of 
the Coulomb matrix elements but also on the frequency 
dependence of U. For TMs with the bcc structure, the 
effective interaction between d electrons in Cg orbitals 
are about 0.3 eV larger than t2g ones, whereas this dif- 
ference is about 0.1 eV for TMs having closed packed fee 
and hep (results not shown) structures. For the former 
TMs, the U{uj) show strong variations at low frequen- 
cies (see Fig. [3]), which suggests that the use of the static 
value U{uj = 0) in model Hamiltonians may be inappro- 
priate. For the latter elements, U (ui) shows a smoother 
behavior. At the plasmon frequency (20-30 eV) the U{Ld) 
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FIG. 3: (Color onhne) (a) U{uj) for bcc V, Nb, and Ta. The 
inset shows J(tj) for the same elements, (b) The same as (a) 
for fee Ni, Pd, and Pt. 
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FIG. 4: (Color online) (a) The ratio U /W of the eflfective 
Coulomb interaction U and the d bandwidth W for the TM 
series. The inset shows the d bandwidths. (b) Density of 
states N{tF) at the Fermi level for the NM states. 



increases rapidly. Above this frequency the screening is 
not effective and U{lo) approaches the bare value (com- 
pare Fig.[l]). In contrast to C/(w), the exchange J{uj) is 
only weakly energy dependent and does not shovif signif- 
icant variations at the plasmon frequency. 

Finally, we discuss the strength of the electronic cor- 
relations and the instability of the paramagnetic state 
towards ferromagnetism for the late Zd TMs. In Fig.|3] 
(a) we show U/W ratios for the TM series, where the 
d bandwidths W are obtained from the single-particle 
band structure. As seen, similar to the bare Coulomb 
interaction V , the U/W ratio increases from early to late 
TMs and, as a result, the correlation strength increases. 
Despite the similar U parameters for isovalent TMs, the 
3d elements have larger U/W ratios than 4d and hd ele- 
ments due to the much smaller bandwidths. For all TMs, 
except Cu and Ag, we have U/W < 1, which reveals 
weak electronic correlations in these materials. Ferro- 
magnetism of the late TMs can be related to large U/W 
ratios, but this condition is not sufficient. In addition, the 
dimensionality and the crystal structure, which dictates 
the shape of the density of states (DOS), is crucial for 
the appearance of itinerant ferromagnetism.^^ In a mean- 



field treatment of itinerant ferromagnetism the instability 
of the paramagnetic state is given by the Stoner criterion 
IN{eF) > 1, where / is the Stoner parameter and N^ep) 
is the DOS at the Fermi level in the NM state. Using 
the Hartree-Fock solution of the multi-orbital Hubbard 
model, StoUhoff et alM proposed a relationship between 
the Stoner parameter / and the Hubbard U and J, which 
is given hy I = {U + 6J)/5. These authors showed that 
electron correlation reduces / by roughly 40%. Using the 
calculated U and J values we get / = 0.98, 1.08, and 1.04 
for Fe, Co, and Ni, respectively, which is very close to the 
values 0.92 (Fe), 0.98 (Co), and 1.02 (Ni) obtained from 
linear-response calculations i^'* Among the 3d series, only 
Mn, Fe, Co, and Ni satisfy the Stoner criterion due to 
the large DOS at the Fermi level [see Fig.|4](b)], and the 
paramagnetic state is unstable towards the formation of 
ferromagnetism. The 4d element Pd is nearly ferromag- 
netic. It shows strong spin fluctuations and exchange 
enhancement . 

In conclusion, by employing a new parameter-free 
cRPA scheme we have calculated the effective on-site 
Coulomb interaction (Hubbard U) between localized d 
electrons in TMs. We have shown that the Hubbard 
U depends on the crystal structure, spin polarization, 
d electron number, and d orbital filling, while it is in- 
sensitive to the d character of the elements. Most of the 
isovalent TMs assume similar U values. The obtained 
U parameters for the 3d TMs are in good agreement 
with previous studies as well as available experimental 
data and predict correctly the paramagnetic instability 
towards the ferromagnetic state for the late 3ds. The 
U (J) values as calculated in the presented approach in- 
crease considerably the predictive power of the LDA-I-J7 
and LDA-I-DMFT schemes applied to describe correlated 
electron materials. 
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